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ABSTRACT 

We analyze multicolor light curves and high resolution optical spectroscopy of the eclipsing 
binai-y and Double Periodic Vai-iable OGLE 05155332-6925581. According to Mennickent et 
al., this system shows a significant change in the long non-orbital photometric cycle, a loop 
in the color-magnitude diagram during this cycle and discrete spectral absorption components 
that were interpreted as evidence of systemic mass loss. We find that the best fit to the multi- 
band light curves requires a circumprimary optically thick disc with a radius about twice the 
radius of the more massive star. The spectroscopy indicates a mass ratio of 0.21 ± 0.02 and 
masses for the hot and cool stars of 9.1 ± 0.5 and 1.9 + 0.2 Mq, respectively. A comparison 
with synthetic binary-star evolutionary models indicates that the system has an age of 4.76 
X 10^ years, is in the phase of rapid mass transfer, the second one in the life of this binary, 
in a Case-B mass-exchange stage. Donor-subtracted Ha profiles show the presence of double 
emission formed probably in an optically thin circumstellar medium, while the variable He I 
profile and the H/3 absorption wings are probably formed in the optically thick circumprimary 
disc. The model that best fit the observations shows the system with a relatively large mass 
transfer rate of M = 3.1 x 10"^ Mo/yr. However, the orbital period remains relatively stable 
during almost 15 years. This observation suggests that the hot-spot mass-loss model proposed 
by other authors is not adequate in this case, and that some other mechanism is efficiently 
removing angular momentum from the binary. Furthermore, our observations suggest that 
the DPV phenomenon could have an important effect in the balance of mass and angular 
momentum in the system. 

Key words: stars: early-type-stars: evolution- stars: mass-loss, -stars: emission-line, bina- 
ries:eclipsing. stars: variables:others 



1 INTRODUCTION 



The investigation carried out by M08 was based on multi- 
wavelength light curve analysis, and resulted in an intermediate- 
mass semi-detached Algol-type binary in an evolutionary stage 
characterized by mass exchange and mass loss where the less mas- 
sive star transfers matter onto the more massive star. M08 also de- 
tected evidence for a relatively luminous disc around the primary 
star. They reported a loop in the color-magnitude diagram during 
the long cycle that was interpreted as evidence of mass loss modu- 
lated by some still unknown mechanism. This was also supported 
by the motion of discrete absorption components observed in the in- 
frared hydrogen lines. M08 offered the explanation of cyclic mass 
loss for the long photometric cycle, but failed in identifying the 
motion of the gainer in their low resolution spectra, and the stellar 




cally the eclipsing LMC star OGLE 05155332-6925581(MACHO 
IDs 79.5739.807 and 78.5739.78; OGLE -SC8-125836). This star 
is a member of the group of close binaries named Double Peri- 
odic Variables (DPVs; Mennickent et al. 2003, 2005a). These stars 
show two closely related photometric periodicities; one longer and 
not exactly cyclic (Mennickent et al. 2005b) and one shorter that 
reflects the binary period. 



* Based on observations carried out at ESO telescope: ESO proposal 

082.D-0575(A). 

t E-mail: hgarrido@eso.org 
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Table 1. Summary of the wavelength coverage, exposure time and spectral 
resolution of the different FLAMES modes and settings used in the obser- 
vation of OGLE 05155332-6925581. R is the spectral resolving power and 
N the number of spectra obtained. 



Mode 


Setting 


/l-range (A) 


Exp. Time (s) 


R 


N. 


Medusa 


HR5A 


4341-4585 


2770 


20000 


14 


UVES 


580 


4775-6800 


2770 


47000 


5 



parameters, in particular the mass ratio, remained relatively uncer- 
tain in their study. 

Consequently, we decided to investigate this system with high 
resolution spectroscopy, in order to look for the missing features of 
the hot star, find definitive stellar parameters and understand better 
the evolutionary stage of the system when comparing with results 
of non-conservative evolutionary models for binary stars. In this 
paper we present our study of this high-resolution spectroscopy. In 
the next section we give a detailed report of the observations used 
for the analysis; our results are presented in Section 3; a detailed 
discussion of these results is given in Section 4 and our conclusions 
are presented in Section 5. 



2 OBSERVATIONS AND DATA REDUCTION 

We have monitored OGLE 05155332-6925581 spectroscopically 
in the optical wavelength range. The observations of this object 
were obtained in the period from 17 January to 26 September of 
2009, using the Fibre Large Array Multi-Element Spectrograph 
(Flames) on the Very large Telescope (VLT) in service mode. In 
total 19 spectra were obtained, 14 of them with the Giraffe spec- 
trograph using the standard Medusa-Giraffe setting HR5A in the 
wavelength region of 4 341 - 4 585 A and 5 spectra using the UVES 
spectrograph in the region 4 775 - 6 800 A with standard UVES 
setting 580, as shown in Table 1. The UVES data do not include 
the region 5 756 to 5 823 A because of the gap between the two 
detectors. 

The ESQ Common Pipeline Library (CPL) FLAMES reduc- 
tion routines were used for initial data processing and consisted 
of flat-fielding, bias subtraction and wavelength calibration. Sepa- 
rate fibers were used to observe the sky in each exposure, each sky 
fibre was inspected for signs of cross-contamination from bright 
spectra/emission lines from adjacent fibers on the detector. Any 
contaminated sky fiber was rejected before creating a median sky 
spectrum, which was then subtracted from each science image. 
Related to the sky subtraction, one of the principal limitations of 
the fiber spectroscopy is the subtraction of local nebular emission. 
Due to the distance of the Magellanic Clouds, even long-slit spec- 
troscopy can suffer difficulties from spatially- varying nebular emis- 
sion jWalbom et al. 2002^ . Because we selected Balmer and He I 
lines for our analysis, we examined these lines for residuals of neb- 
ulosity subtraction. These were actually found near the Ha profile, 
as shown later, but the affected regions are easily identified and iso- 
lated, and not considered in our analysis. The above method was 
implemented by [Evans et al. (2011)| to reduce VLT-Flames obser- 
vations. 

The UVES frames are reduced using a similar method as for 
the Giraffe data. We removed cosmic rays and then extracted the 
spectra, that were normalized over the entire spectral region. All 
spectra discussed in this paper are corrected to the heliocentric 
frame, using the RVCORRECT and DOPCOR packages imple- 



Table 3. Principal lines detected in the optical spectrum of 
OGLE 05155332-6925581. The symbol t indicate lines used in the 
RV measured. 



^(A) 


Line 


A{k) 


Line 


A{k) 


Line 


4 387.9 


He I 


4471.5t 


He I 


4 390.6 


Mg 11 


4481.2t 


Mgll 


4 549.5 t 


Fell 


4861.3 


H/3 


4921.9t 


He I 


5 015.7 


He I 


5 875.6 


Hel 


5889.9 


Na I int. 


5 895.2 


Na I int. 


6562.8 




6678.2 


He I 











mented in IRAlQ We use indifferently the words primary {gainer) 
and secondary {donor) for the components of the binary system as 
usual in the literature of mass exchanging semi-detached binaries. 



3 RESULTS 

3.1 Radial velocity measurements and spectroscopic 
mass-ratio 

M08 could not estimate a spectroscopic mass ratio for this system, 
mainly due to the weakness of the features of the primary star in 
the available low resolution spectra. In addition, the analysis of the 
OGLE II and MACHO light curves provided inconsistent values 
for the photometric mass ratio, in part since the light curve model 
at that time did not include a contribution from a potential circum- 
primary disc. In this section we provide the first detection of the 
features of the hotter component and an improved determination 
for the the system mass ratio. 

Line identification was done with the aid of the spectral atlas 
for B main-sequence stars provided on-line by Heidelberg Univer- 
sit)]^ A list of the most prominent lines identified in the spectra 
are given in Table 3. We identified several strong absorption lines 
of He I, Mg II and Fe II in our Flames spectra. The UVES spectra 
display mainly Balmer and strong He I lines. Ho- profiles are broad 
absorptions filled by emission. In this line we detected contribu- 
tions of both stellar components plus a very narrow central absorp- 
tion spike. This feature turned to be residual absorption from the 
process of sky subtraction, consistent with spatially variable neb- 
ular emission in this part of the sky. The He I lines of the primary 
appear to be partially blended with lines from the secondary star, 
indicating a relatively hot donor star. 

All spectra were normalized to the continuum before perform- 
ing the RV measurements and analysis. The RVs were measured by 
calculating the positions of the line centre with a gaussian fit using 
the IRAF SPLOT task. We measured the heliocentric radial veloc- 
ities for the secondary using the Mg II 4 481 A and Fell 4 549 A 
lines. These RVs are listed in Table[2] along with the corresponding 
orbital phase computed with the ephemerides given in Section [T2| 
and the errors taken as the standard deviations of the measurements. 
All velocities are considered with equal statistical weight. We as- 
sumed a circular orbit, so a sine function was selected to fit them, 
obtaining a systemic velocity ydonm = 253.20 ± 1.22 kms"' and 



IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in 
Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 

^ http://www.lsw.uni-heidelberg.de/cgi-bin/websynspec.cgi 
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Table 2. RVs measured from the Medusa and UVES spectra for OGLE 05155332-6925581referred to the Local Standard of Rest. The orbital phase <Do is 
calculated following the ephemeris given in Eq. (1). Subindex 1 and 2 refers to the gainer and donor, respectively. 



HJD 


<f>o 


RVi 


(O-C)i 


RV2 


(0-C)2 


Mode 


-2450000 




km s"' 


km s"' 


km s"' 


km s"' 




Jl./O'T. OO 7 J 


043 


z.^o.oz. ni Z..UJ 


_19 ^9 + 9 6^ 


'^OS 71 + 118^ 

JLf J. / I X 1 . 1 O 


—0 61 + 1 18 




S 1 no R7Q1 


W.Z JO 


1 o / .to i i .UO 


—8 44 -(- 1 68 


/ISS 9Q -(- 1 ^7 


2 28 + 1 37 


^dedusa 


to J J .uuo / 


292 


178 + 70^^ 
1 1 o.jy X w. / 7 


— 11 21 + 79 






UVES 


to J J .OJ J o 




1 87 Q4 + 1 '\C\ 


— 1 45 + 1 50 


444 18 -t- n 7S ^ 

^Ht. JO i U. / J 


—0 14 + n 7S 


^dedusa 


t O U O . VJ\7 J VJ 


351 


1 Q7 64 + n 46 


Q + n 46 






UVES 


5079.8642 


0.353 


185.83 ±0.94 


-2.55 ± 0.94 


407.31 ±0.48* 


-5.67 ±0.48 


Medusa 


4868.6495 


0.357 


183.77 ± 1.24 


-4.73 ± 1.24 


417.92 ± 1.50* 


8.04 ± 1.50 


Medusa 


4876.5457 


0.441 


215. 14± 1.89° 


18.33 ± 1.89 






UVES 


4898.6029 


0.469 


207.32 ±0.51 " 


5.56 ±0.51 






UVES 


5080.8272 


0.485 


215.98 ± 1.54 


10.98 ± 1.54 


265.49 ± 1.02* 


-6.14 ± 1.02 


Medusa 


4848.5946 


0.604 


228.10 ±0.98 


-6.29 ± 0.98 


133.30 ±0.34 


2.69 ± 0.34 


Medusa 


4848.6366 


0.610 


245.81 ± 1.25" 


9.87 ± 1.25 






UVES 


5081.8340 


0.623 


223.96 ± 1.41 


-15.33 ± 1.41 


116.96 ±0.72 


4.48 ± 0.72 


Medusa 


4892.5550 


0.639 


233.16 ±0.98 


-10.16 ±0.98 


100.81 ± 1.21* 


2.06 ± 1.21 


Medusa 


4863.5947 


0.663 


238.02 ±3.16 


-11.07± 3.16 


80.23 ± 1.33 


-0.88 ± 1.33 


Medusa 


4907.5324 


0.695 


261.07 ± 1.81 


4.98 ± 1.81 


61.32 ±0.16* 


-2.39 ±0.16 


Medusa 


5082.8394 


0.761 


264.24 ± 1.80 


-2.68 ± 1.80 


56.16 ±0.60* 


3.85 ± 0.60 


Medusa 


4865.5906 


0.937 


280.54 ±2.76 


17.01 ±2.76 


180.51 ±0.88 


5.29 ± 0.88 


Medusa 


5076.9004 


0.946 


272.10 ± 1.73 


10.02 ± 1.73 


180.56 ± 1.58* 


-5.26 ± 1.58 


Medusa 



Note. " RVs values based on the line He II 4 921 A, * RVs values based on the lines Mg II 4481 A, only. 

semi-amplitude A^donoi = 200.42 ± 1.66 kms"' with rms 4 kms"' 
(see Fig.[T](. 

The lines He 1 4 47 1 A and He 1 4 92 1 A mainly follow the pri- 
mary star, but a minor contribution from the secondary star is 
present. We used a deblending method to isolate the components 
in the analysis. We performed this task with the IRAF deblending 
'd' routine available in the SPLOT package. We carefully measured 
the RVs of the He I lines (see Table|2j and fit them with a sine func- 
tion, obtaining a systemic velocity yjaincr = 229.85 ± 2.87 kms"' 
and a semi-amplitude K^.^[^^„ = -41.69 ± 3.67 kms"' with rms 
11 kms"'. Although there are other He I and Balmer lines in our 
UVES spectra with similar double profiles around phases 0.29 and 
0.61, we could not de-blend them to include them in the radial ve- 
locity calculation, since the individual components were hardly re- 
solved. The difference in the systemic velocity between the solu- 
tion for the primary and secondary can be explained by emission of 
circumstellar matter affecting the He I lines. From the above semi- 
amplitudes we find a spectroscopic mass ratio of q^p = 0.21 + 0.02 
for OGLE 05155332-6925581, which is smaller than q = 0.29 
obtained by M08 from the best MACHO blue light-curve fitting 
model. 

Interestingly, our donor half-amplitudes derived from the 
study of Fe II and Mg II lines are larger than those found by M08 
in the Balmer lines by up to 80 km s ' . This could be associated to 
a larger contamination by emission in Balmer lines, as suggested 
by our donor-subtracted Ha profiles (Section 3.3). The fact that the 
emission formed around the gainer moves in anti-phase with the 
donor absorption produces the net effect of lowering the observed 
donor radial velocity half-amplitude. We considered another possi- 
ble cause for the half-amplitude discrepancy; displacement of the 
light-centre due to donor irradiation by the disc and gainer (e.g. 
Wade & Home 1988, Watson et al. 2003). This efi'ect would re- 
sult in an over-estimate of K2, leading to a corresponding under- 
estimate of q. However, the absence of asymmetries in the donor 
line profiles and the large observed differences in HI line radial 




Figure 1. The RV curve together with O-C for OGLE 05155332-6925581. 
Plot is folded to orbital period using the ephemeris for the main minimum, 
given in the equation^ Panel shows the RV measurements with the best-fit 
solutions. Filled circles represent the gainer and triangles the donor. The 
bottom panel shows the residuals. 



velocity half-amplitudes (up to 45 kms"'), suggest that residual 
emission is the most likely explanation. It is difficult to imagine 
how light-centres could difi"er so much between lines of the same 
element at a given ionization stage. 
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Figure 2. Observed (LCO) and synthetic (LCC) light-curves of LMC star OGLE 05155332-6925581 obtained by analyzing MACHO blue, red and OGLE 
Lband photometiic observations; final O-C residuals between the observed and optimum synthetic light curves; fluxes of donor, gainer and of the accretion 
disc, normalized to the donor flux at phase 0.25; the views of the optimal model at orbital phases 0.25 and 0.75, obtained with parameters estimated by the 
light curve analysis. 



3.2 The light curve analysis and model of the accretion disc 

We have used public domain photometric data available from the 
MACHO ( |Alcock et al., 1993| l and OGLE II and III databases 
( |Udalski et al., I997||Szymanski 2005| and |Poleski et al. 2010^ . The 



photometric data cover 13.7 years. The observations were folded to 
the orbital period using the ephemeris for the main minimum from 
M08: 



Tn,in(HJD) = 2450000.1392(21) -t- 7''.284297 (10) x E. 



(1) 
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Table 4. Result of the analysis of LMC star OGLE 05155332 - 6925581 in MACHO blue, MACHO red and OGLE I-band light curves obtained by solving 
the inverse problem for the Roche model with an accretion disc around the more-massive (gainer) component in critical rotation regime. 



Quantity 


MACHO blue 


MACHO red 


OGLE I 


mean BRl 


Quantity 






2964 


2 796 


883 




Ml r A4^1 


9 1 + 05 




1 .9490 


1.8550 


0.2634 






1.9 ± 0.2 




0.0256 


0258 


0173 




'^/jtJ^Oj 


5 6 + 02 


'L J 


82.33 


82.25 


82.36 


82.3 ± 0.3 


K IR^l 


8.9 ± 0.3 




76 


77 


77 


77 + 0^ 


log 


3 9 + 01 


T ,rKi 


12 200 


12 800 


12 770 


1 9 600 + 600 


log Sc 


2 8 + 01 


*^t"orbJ 


020 


021 


025 


09 -(- n 1 


™bol 


—53 + 02 


dc[«orb] 


0.070 


0.072 


0.077 


0.07 ± 0.01 




-3.4 ±0.1 


aT 


4.0 


3.5 


3.5 


3.7 ± 0.3 


Oorb[Ro] 


35.2 ±0.5 


fh 


10.6 


10.6 


10.5 


10.6 ±0.3 


«rf[Ro] 


14.1 ±0.5 


Fh 


1.00 


1.00 


1.00 


1.00 


de[Ro] 


0.8 ± 0.3 


Tc[K] 


12 930 


12840 


12 830 


12 900 ±500 


dc[Ro] 


2.6 ±0.5 


Ahs =Ths/Trf 


1.13 


1.12 


1.10 


1.12±0.1 








14.1 


16.3 


16.2 


15.5 ± 1.0 








327.4 


322.9 


322.4 


324.2 ±5.0 






erad[°] 


-34.2 


-33.8 


-34.9 


-34.3 ± 5.0 






Ate =Tte/Trf 


1.15 


1.19 


1.20 


1.18±0.1 






ebs[°] 


54.1 


52.2 


48.2 


51.5 ±3.0 






^bs[°] 


120.5 


136.5 


134.2 


130.4 ± 13.0 








7.94 


7.92 


7.90 


7.92 ± 0.02 






He 


2.56 


2.56 


2.56 


2.56 ± 0.02 







FIXED PARAMETER: q = Md Mi, = 0.21 - mass ratio of the components, Ti, = 25 000 K - temperature of the more-massive (gainer), F^ = 1.0 - filling 
factor for the critical Roche lobe of the donor, Fh = R;,/R;r = 1.0 -filing factor for the critical non-synchronous lobe of the more massive gainer (ratio 
of the stellar polar radius to the critical non-synchronous lobe radius along z-axis for a star in critical rotation regime), = 1.00 -synchronous rotation 
coefficients of the donor, /Jh.c = 0.25 - gravity-darkening coefficients of the components, Ah,c = 1.0 - albedo coefficients of the components. 
Note: n - total number of observations in MACHO blue (B), red (R) and OGLE I-band; S(0 - Cf - final sum of square of residual between observed 
(LCO) and synthetic (LCC) light -curves, (Ti-ms - root-mean-square of the residuals, / - orbit inclination (in arc degrees), F^i — Rj/Ryj; - disc dimension 
factor (the ratio of the disc radius to the critical Roche lobe radius along y-axis), T;j - disc-edge temperature , de, dc - disc thicknesses (at the edge and at 
the centre of the disc, respectively) in the units of the distance between the components, ax -disc temperature distribution coelficient, fh - non-synchronous 
rotation coefficient of the more massive gainer (in the critical rotation regime), Tc -temperature of the less massive cooler donor, flh.c - dimensionless 
surface potentials of the hotter gainer and cooler donor, A1/,_r[A1o], 'K;, ^[Ro] - stellar masses and mean radii of star in solar units, loggj ^ - logarithm (base 
10) of the system components effective gravity, M^J^'j - absolute stellar bolometric magnitudes, aorb[Ro]> '^(/[Rq]. de[RG]. dc[R<3] - orbital semi-major axis, 
disc radius and disc thicknesses at its edge and centre, respectively, given in solar units. 



The light-curve analysis was performed using the inverse-problem 
solving method l |Djurasevic 1992) based on the simplex algorithm 
and the current version of the model of a binary system with a 
circumprimary disc ( [Djurasevic et al. 2008) . The model and code 
have been widely used and tested during our recent research of 
intermediate-mass interacting binaries (e.g. Djurasevic et al. 2010, 
2011, 2012, Mennickent et al. 2012a). For our calculations we as- 
sume that the disc is optically thick and its edge can be approxi- 
mated by a cylindrical surface. The thickness of the disc can change 
linearly with radial distance, allowing the disc to take a conical 
shape (that it can be convex, concave or plane-parallel). The geo- 
metrical properties of the disc are determined by its radius (7?^), its 
thickness at the edge (de) and the thickness at the centre (dc). 

The cylindrical edge of the disc is characterized by its tem- 
perature, Tj, and the conical surface of the disc by a radial tem- 
perature profile obtained by modifying the simple a temperature 
distribution. Assuming that the disc is in physical and thermal con- 
tact with the gainer, so the inner radius and temperature of the disc 
are equal to the temperature and radius of the star {%„ T;,). In this 
way, we describe the temperature distribution by the following law 
( [Djurasevic et al. 2008| i: 

nr)-Td^[n-T.)[^-[^J\ (2) 



The temperature of the disc at the edge (T^) and the temperature ex- 
ponent (ax), as well as the radii of the star (7?/,) and of the disc (7?^) 
are free parameters, determined by solving the inverse problem. 

According to the temperatures of the system components, 
we set the gravity-darkening exponents and component's albedos 
at their theoretical values (j8i 2 = 0.25, Aj 2 = 1.0), corresponding 
to von Zeipel's law for fully radiating shells and complete re- 
radiation. 

Limb-darkening was applied by using the nonlinear approx- 
imation by |Claret| ( |2000| l for the photometric bands. The limb- 
darkening coefficients are interpolated in each iteration for the cur- 
rent values of the effective temperature, Tefj, and surface grav- 
ity, log[Q^, as is described in detail by [Djurasevic et al. (2004)| 
The limb-darkening was applied to the disk in the same way, with 
logjo g that corresponds to the middle of the disk radius. 

The models of the system can be refined by introducing active 
regions on the edge of the disc. The active regions have a higher 
local temperatures so their inclusion results in a non-uniform dis- 
tribution of radiation. The existence of such regions (hot and bright 
spots) can be explained by the gas dynamic of interacting binaries, 
see e.g. [Heemskerk[ ( |T994l l and Bisikalo et al. (1998, 2000, 2005). 
Based on these investigations of systems with accretion discs, we 
recently updated our model to allow inclusion of the hotspot (hs) 
and bright spot (bs) active regions. 
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These regions are characterized by the following parameters: 
Ahs.bs =Ths.bs / Trf - the hot and bright spot to disc temperature ra- 
tios, Shs.bs and /Ihs.bs - the spot's angular dimensions and longitudes 
(in arc degrees) and Sj^d - the angle between the line perpendicu- 
lar to the local disc edge surface and the direction of the hot-spot 
maximum radiation. The longitude A is measured clockwise (as 
viewed from the direction of the +Z-axis, which is orthogonal to 
the orbital plane) with respect to the line connecting the star cen- 
tres (-nX-axis), in the range - 360 degrees (for more details, see 
[Djurasevic et al., 2010^ . These parameters are also determined by 
solving the inverse problem. Including these active regions in the 
model, we found significant improvements in the quality of the fit 
to observations. For further information about these active regions 
in the model, see [Djurasevic et ar~| ( |2012| . 

One possible limitation of the code at its present implementa- 
tion is the lack of a detailed treatment for the donor irradiation by 
the disk part facing the donor, including hot spot. This effect could 
be potentially important in close binaries with a large difference 
in stellar and disk temperatures. However, it should be a second- 
order effect compared with the stellar and disc flux contributions 
(including donor irradiation by the gainer) already implemented in 
our code. This is demonstrated in the very good fit obtained with 
the orbital multi-wavelength light curves. 

3.3 Fitting Procedure for photometric data 

The fitting was done for a limited set of important parameters. The 
mass ratio q = MdMi, =0.21 (where the subscripts h and c stand 
for the more massive hotter gainer and less massive cooler donor) 
was set to the value determined by the radial velocities of our new 
FLAMES observations, discussed in Section [TT] The temperature 
for the primary was fixed to the value given by M08, viz. 25 000 
K. Rotation for the donor was assumed synchronous (fc = 1.0) 
since it is assumed that it has filled its Roche lobe (i.e. the filling 
factor of the donor was set to = 1.0). In the case of the gainer, 
however, the accreted material from the disc is expected to transfer 
enough angular moment to increase the spin rate of the gainer up 
to the critical rotation velocity as soon as even a small fraction of 
the mass has been transferred ( |de Mink , Pols & Glebbeek 2007[ l. 
This means that the gainer fills its corresponding non-synchronous 
Roche lobe for the star rotating in the critical regime. This assump- 
tion is justified by the rapid spin-up of the gainer at the system evo- 
lutionary stage (see Section 4). So, the dimensions and the amount 
of rotational distortion are uniquely determined by the factor of 
non-synchronous rotation, which is the ratio between the actual 
and the Keplerian angular velocity. In this model the factor of non- 
synchronous rotation for the gainer was a free parameter also deter- 
mined by solving the inverse problem. However, we also calculated 
the best fit model for a synchronous gainer, and the results practi- 
cally do not differ from the critical case. The results of our analysis 
are presented in Table 4. 

The stellar parameters improve significantly the solution given 
by M08. We find primary and donor stars with 9.1 ± 0.5 Mq and 
1.9 ± 0.2 Mq, respectively. They have log g of 3.9 ± 0.1 and 2.8 
± 0.1 and radii of 5.6 ± 0.2 Rq and 8.9 ± 0.3 R^, respectively. 
The bolometric magnitudes are M'^^ = -5.3 and M[J^j = -3.4. 
The best fit requires an optically thick disc around the gainer. This 
disc has a moderately convex shape, with central thickness d;. a; 
2.6 Rq and thickness at the edge de « 0.8 Rq. The radius of the disc 
is X 14^o. more than twice larger than the gainer radius and 
about ~ 77% of the Roche lobe radius. The temperature of the disc 
increases from T^; = 12 600 K at its edge to T;, = 25 000 K in the 



inner radius, where it is in thermal and physical contact with the 
gainer. 

The two spots on the edge of the disc mentioned in Section 
3.2 were found at longitudes 324 and 130 degrees with tempera- 
tures 12 and 18 % higher than the disc edge temperature, respec- 
tively. Whereas the hot spot represents the region of impact of the 
gas stream with the disc, the bright spot atAx 130 degrees can be 
interpreted as a region where the disc significantly deviates from 
circular shape, or as a region where the stream material falls back 
to the opposite side of the disc, after being deflected by the impact. 
Due to the relatively small size and temperature difference com- 
pared to the lateral side of the disk, the hot spot is not sufficient to 
irradiate the face of the donor significantly. The bright spot is lo- 
cated on the disk lateral side which is not able to contribute to the 
irradiation of the donor. 

In Fig.|2]we show the comparison between the synthetic and 
the OGLE-I, MACHO-B and MACHO-R light curves. We do not 
include errors bars since the scatter in the light curves probably 
gives a better representation of the true errors. In general the fit is 
pretty good, but there is a slightly larger scatter in O - C around 
secondary eclipse (I-band) and main eclipse (b and r bands). This 
scatter is very minor and could represent a third light in the sys- 
tem that is not considered in the model and that is best visible 
during eclipses. A probable explanation for these larger residuals 
is the 3D structure of the hot spot, including material located at 
high latitudes, whose small contribution to the total flux is detected 
only during eclipses. We show also the geometrical model of the 
system at quadratures. The residuals do not depend on the long- 
cycle phase. This important result indicates that the orbital solution 
including the disc remains constant during the whole long cycle. 
Therefore, the disc represented in the model is not the cause for the 
long cycle. The figure also illustrates the geometrically thin nature 
of the disc, its relative size and the positions of the hot and bright 
spots in the system. At quadratures, the disc contributes about 15% 
to the total light of the system at the /-band, whereas at main eclipse 
the secondary contributes about 84%. 

3.4 Donor-subtracted and residual spectra 

One of the goals of this paper is to find information about the en- 
vironment in the vicinity of the hot primary star. However, from 
the disentangling processes we know that our composite spectrum 
is contaminated by the light contribution of the B-type donor, in 
particular in the Balmer and He lines. Thus, we need to remove the 
secondary star's spectral contribution from the observed spectra to 
study the light contribution of the primary and its circumstellar gas. 

In order to remove the donor light from the composite spec- 
tra we assume that its contribution to the total light is additive to 
the other light sources and represented by the light-curve model 
proposed in the section [X2] Furthermore, we compute a synthetic 
spectrum using the data in Table 4 for the donor star. 

The model for the donor stellar atmosphere was constructed 
for an effective temperature of 13 000 K and surface gravity of 
logg = 3.0 with solar abundance, using the LTE ATLAS9 code 
l |Kurucz, 1993^ , which handles the line opacity with the opacity 
distribution functions method (ODFs). The Kurucz models are con- 
structed with the assumptions of plane-parallel geometry and hy- 
drostatic and radiative equilibrium of the gas. The synthetic spec- 
trum was computed with the SYNTHE code ^Kurucz, 1993[ l and 
broadening at a rotational velocity of Vrot sin ; = 60 km s ' , as- 
suming that the donor has synchronous rotation, as was mentioned 
in Section 3.3. Both codes, ATLAS9 and SYNTHE were ported 
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Figure 3. Donor-subtracted spectra showing Yip variable extended wings up to ±1000 km s"'. He I 4921 A absorption line is also visible. The red line 
overplotted on the donor-subtracted spectra is the SYNTHE model using Kurucz atmospheres for a more-massive (gainer) component in critical rotation 
regime. The labels in the upper right of each spectrum indicate the orbital phase. 
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Figure 4. Residual spectra showing the Ha and He I 6 678 line profiles after removing the donor and gainer synthetic spectra. The nan'ow absorption in Ha is 
an artifact of the reduction process (Section 2). The labels in the upper right of each spectrum indicate the orbital phase. 



tmder GNU Linux by |Sbordone et al.,H2005^ and are avail able on- 
lin^ The atomic data were taken from |Castelli & Hu brig ( 2004p] 
After that, the synthetic spectrum was Doppler corrected and 



3 wwwuser.oat.ts.astro.it/atmos/ 

wwwuser.oat.ts.astro.il/castelli/grids.html 



scaled according to the contribution of the donor star at a given 
orbital phase and at the given wavelength range. This method was 
used by [Mennickent et al.] j2012b( l to remove the donor spectrum 
from the observed composite spectrum of the interacting binary 
V393 Sco. Then the isolated gainer spectrum plus disc in units of 
its normalized flux is given by 
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fa (/I, Do, = m Do, - P (Do, 4) x /d {A, Do) 



(3) 



where /g is the donor-subtracted flux, / the observed flux, /d the 
synthetic donor spectrum, P the fractional contribution of the donor 
derived from our hght-curve model and the representative wave- 
length where P is calculated. The theoretical light contribution fac- 
tors for the donor at different spectral lines according to the light- 
curve model are shown in Fig|5] The result of this process was a 
set of "donor-subtracted" spectra that were normalized to the new 
continuum. 

Our spectra cover 1.5 times the long cycle of 172 days, but 
only the Medusa spectra are well distributed to adequately map the 
long variability. We inspected these profiles, having subtracted the 
donor, but without observing significant changes in the small wave- 
length range provided by the Medusa spectra. In the following, we 
analyze the UVES spectra, that show the prominent Hy6 and He I 
4 921 A lines, looking for orbital variability. 

The donor-subtracted spectra for the H/3 and He 1 4 921 A ab- 
sorption lines in orbital phases 0.29 to 0.61 are shown in Fig.|3] 
The Hy6 profiles are broad and show variable absorption wings ex- 
tending up to ±1000 km s"'. These lines are more symmetric near 
the secondary eclipse. The Helium lines are broad and asymmet- 
ric and sometimes show flanking emission, especially at the red 
side of the He 1 4 921 line in Fig. 3. Measures of Full Width at Half 
Maximum (FWHM) and equivalent width (EW) for these lines are 
shown in Table 5 for each orbital phase. Fig. 3 also shows a syn- 
thetic spectrum computed with SYNTHE for the gainer using pa- 
rameters Teff = 25 000 K and logg = 3.9 and solar abundance. Our 
synthesized spectrum was computed for a gainer in critical rotation 
using a rotational parameter V,,,, sin i = 412 km s"' , as estimated in 
Section l33] 

It is clear that the profiles are highly variable in shape. The fact 
that the H/3 line is deeper than the gainer synthetic profile suggests 
that it is produced in an extended absorption media beyond the stel- 
lar photosphere, a kind of pseudo-photosphere. This media could 
be the disc revealed in the light curve analysis. The lack of obvi- 
ous variability in the 13.7 year folded orbital light curve indicates 
a stable disc, so the large profile variability suggests the presence 
of additional mass streams in the system. The general features are 
compatible with absorption/emission in an extended photosphere 
around the gainer plus additional mass flows. 

We found a good match using the critical rotational velocity 
for the He I profiles and Hj8 absorption wings of our subtracted 
spectra, around secondary eclipse, when the gainer is in front to the 
donor, as shown in Fig. [3] However, we observe a small deviation at 
phase 0.29, when the wings of the line profiles are less symmetrical, 
showing extended absorption in the H/3 red wing and emission in 
the He 4 921 A red wing. 

We obtained residual spectra around the Ha region. These 
spectra were obtained from the donor-subtracted spectra after re- 
moving the gainer synthetic spectrum. The process reveals weak 
double emission at Ha with a peak separation of about 620 km s"' 
(Fig. 4). The same exercise at the He I 6 678 A line shows absorp- 
tion and traces of emission. We are cautious about interpreting the 
Ha double emission as produced in the circumprimary disc, since 
in the case of V 393 Scorpii the evidence indicates that this kind of 
profile is produced in a bipolar wind (Mennickent et al. 2012b). 

We conclude that line profile fitting cannot be used as proof 
of critical rotation in this system due to the high variability of the 
profiles. On the other hand, the detection of Balmer and He I emis- 
sion is consistent with the presence of circumstellar matter in the 
system. 




Orbital Phase 

Figure 5. Light contribution factor for the donor at different .spectral lines 
and orbital phases according to the light-curve model given in Section [X2] 
L[ , Li and Lj are the gainer, donor and disc fluxes, respectively. 



FWHM (km s"') 



Phase 



H, 



He I [4921] 



EW(10"^A) 
H/3 He I [4921] 



0.292 


579.30 ± 3.60 


335.80 ± 1.08 


444 ± 1 


73 ± 1 


0.351 


574.00 ± 1.19 


338.10 ±0.24 


432 ± 1 


63 ±3 


0.441 


462.53 ± 1.05 


301.73 ±4.15 


364 ± 1 


54 ±2 


0.469 


523.13 ± 1.90 


247.00 ± 3.53 


333 ± 1 


40 ± 1 


0.610 


562.40 ± 0.65 


233.30 ±4.41 


442 ± 1 


58 ± 1 



Table 5. Equivalent width and FWHM of the 1^8 and He 4 92 1 A absorption 
lines from the donor- subtracted UVES spectra. 



4 DISCUSSION 

4.1 On the evolutionary stage of OGLE 05155332-6925581 

The comparison of the stellar parameters of OGLE 05155332- 
6925581 with predictions of evolutionary models for single stars of 
solar metallicity indicates that the gainer is located near the main- 
sequence band, a little displaced to higher luminosities for a star 
of comparable mass, while the donor is evolved from the main se- 
quence, much more luminous for its mass, a fact explained by its 
expansion and subsequent filling of its Roche lobe (see Fig. [6]l. The 
small luminosity excess of the gainer can be explained in terms of 
their slightly advanced evolutionary stage; it has consumed about 
half of the hydrogen in its core, as will be shown later in this sec- 
tion. 

In order to understand the evolutionary stage of this binary we 
compared our system parameters with those of the 561 conservative 
and non-conservative evolutionary tracks given by |Van Rensbergen] 
|et al.| ( |2008l l, available at the Center de Donnees Stellaires (CDS). 
Similarly to the case of V393 Sco studied by Mennickent et al. 
(2012a, hereafter M12a), a multi-parametric fit was carried out for 
each synthetic model (labeled with ;) through the evaluation of the 
qaaniity xj j defined by: 



-y 



(JJk 



S,,.k-Ot 



(4) 



where A' is a normalization factor and oj^- the statistical weight of 
the parameter Ou-, calculated as: 
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Figure 6. Comparison between evolutionary tracks for single stars 
with solar metallicity (Claret 2004) and the physical parameters of 
OGLE 05155332-6925581 in the logTeff- logL diagram. Derived mass 
and logg values are given between parentheses. The evolutionary tracks 
are labelled with the initial masses. 
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Figure 7. Evolutionary tracks for the binary star model from Van Rensber- 
gen et al. (2008) that best fit the data. Donor (red, vertical dashes track) and 
gainer (blue, horizontal dashes track) evolutionary paths are shown, along 
with the observations for OGLE 05155332-6925581(filled circles with er- 
rors bars). The best fit is reached at the time corresponding to the model 
indicated by stars, that is characterized in Table 6. Stellar sizes are propor- 
tional to the dash sizes. 



^Otle(Ot) 



(5) 



where e(Ot) is the error associated with the observable 0^. Sjj^i^ is 
the theoretical parameter for the observable Ok, in the time tj for the 
i-labeled model. For further information about the fitting procedure 
see Ml 2a. 

We find the absolute x' minimum in the tidal-strong interact- 
ing model with initial masses of 8 and 3.2 Mq and initial orbital 
period Po,bj = 2.5 days. The absolute minimum ^^^^ identifies the 
current age of the system along with the theoretical stellar and or- 
bital parameters, which are presented in Table 6. 

The best fit indicates that OGLE 05155332-6925581 has an 
age of 4.76 x 10' years. The corresponding evolutionary tracks for 



the primary and secondary stars are shown in Fig.|7] along with the 
position for the best model for our system. 

We notice a good agreement in the evolutionary tracks given 
by the best model. The secondary star is inflated {R^ = 9.1 Rg) 
and exchanges mass through the LI Lagrangian point at a rate of 
M = 3.105 X lO^^Mo^r"'. This rather large mass transfer rate 
compares relatively well with the low limit found by M08, viz. 8.6 
X 10"' Mg/yr, derived assuming accretion powered luminosity. 



4.2 On mass flows and angular momentum loss 

The orbital light curve was re-examined with the program Pe- 
riodO^ We calculated the error in the main frequency of the fit to 
the orbital light curve. The error consistently given by Monte Carlo 
simulations and the method of least squares is 2 x 10"' Hz. This 
means that the period could drift at most by 2 x 10"' d in 3365 d 
(the OGLE dataset time baseline). This implies that a period change 
at constant rate, if present, should be less than 0.2 s/yr to be con- 
sistent with the photometric time series. However, according to the 
best theoretical model, we find the system in a stage characterized 
by a rapid change of the orbital period (Fig. 8). The period should 
be changing by 4.7 s/yr. This apparent inconsistency was already 
noted by M08, who argued that the effects of mass exchange and 
mass loss could be balanced in the system producing a net effect of 
non-variability for the orbital period. 

The models by Van Rensbergen et al. (2008, 2011) 
parametrize mass and angular momentum loss from the system 
through the parameters /? and r;, respectively. The mass loss is 
driven by radiation pressure from a hot spot located in the stream 
impact region on the stellar surface or accretion disc edge. The 
mass loss extracts angular momentum from the system and in this 
particular model only from the gainer. According to the models of 
Van Rensbergen et al. (201 1), particular pairs of the parameters /3 
and T] should result in a constant orbital period. 
The mass and angular momentum loss from a system during mass 
transfer can be described with the (fi, ;7)-mechanism (see Rappaport 
et al. 1983). Here, 



Mr 



(6) 



is the mass transfer efficiency, i.e. the fraction of mass lost by the 
donor accreted by the gainer. The angular momentum loss J from 
the system resulting from a given amount of mass loss M is deter- 
mined by rj, according to 



j= V^M, (l-yS) 



(7) 



with a the orbital seperation and P the orbital period. It can be 
shown (see e.g. Podsiadlowski et al. 1992) that whenever p and r) 
remain constant throughout the considered timeframe (with indices 
i for initial and f for final), the resulting period variation is given 
by 



forO<j8< 1: 
Pf_ ^ / Mrf-HM,,f 

foryS = 0: 



Mrf 



(8) 



http://www.univie.ac.at/tops/period04/ 
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Table 6. The parameters of the Van Rensbergen et al. (2008) model that 
best fits the OGLE 05155332-6925581data. The hydrogen and helium core 
mass fractions {X^ and Yc) are given for the cool and hot star 
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Age 


4.76 X 10' yr 


Period 


7.241 d 
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1.945 Mo 
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9.249 Mo 




Mc 


-3.17 X IQ-^Meyr-^ 


Ml, 


3.105 X lO-^MoVr-l 




iogTc 


4.120 K 


log Ti, 


4.385 K 




logic 


3.351 Lo 


logi;, 


3.862 Lo 


p 
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4.821 Rq 
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0.00 
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Figure 9. and /? parameters controlling the mass and angular momentum 
loss from the binary according to the model proposed by Rensbergen et al. 
(201 1). The space between the upper and lower curves is for systems with a 
constant orbital period, according to the observed boundaries for the orbital 
period variability. The region above the upper curve is for systems showing 
an orbital period increase, while that below the lower curve is for systems 
with decreasing orbital period. The value for OGLE 05155332-6925581 
of 8.9 X lO"'' is incompatible with a constant orbital period under the view 
of the aforementioned model. 



Pi 



'M,f + Mi,s 



3(V^-l) 



"hi 



(9) 



If it is assumed that mass is lost with the specific orbital angular 
momentum of the gainer, it can be shown that this yields 



Mr 



[ + Ml, 



(10) 



resulting in ;/ <c 1. On the other hand, if one assumes that matter 
is lost by the formation of a non-corotating circumbinary disc af- 
ter passing through the second Lagrangian point, it was shown by 
Soberman et al. (1997) that a typical value is 77 = 2.3. 

Given the observational upper limit for the period variation 
in this paper, we have calculated the maximum amount by which 
the period may have changed during those 15 years. Then, with the 
equations for the period variation during (non-)conservative mass 
transfer and using the physical parameters of this particular system, 
we have calculated that 15 year period variation theoretically for all 
different combinations of yS and j]. Results are shown in Fig|9] 

The zone between the two curves is where the orbital period 
does not change more in 15 years than is allowed by the observa- 
tions. The zone above the upper curve is where the period increases 
too much. This is the case for large /3 (little mass loss) and/or small 
J] (little angular momentum loss). The zone below the lower curve 
is where the period decreases too much. This is the case for small 
p (much mass loss) and large 77 (much angular momentum loss). 

The assumption of specific gainer orbital angular momentum 
loss (the one used in the calculation of the models) yields t] = 
0.00089 for this system. It is obvious that this is way too low in 
order to fall into the constant period region, and that under this as- 
sumption, irrespective of /?, the period will increase (much) more 
than allowed by the observations. Based on this result we argue that 
another source of angular momentum loss is present in the system, 
possibly the outflows through L2 and L3 reported by M08. Notice 
that for r] x 2.3, representative for angular momentum loss from a 



circumbinary disc (Soberman et al. 1997), there is a mass loss (fi w 
0.5) compatible with period constancy in this system. 



5 CONCLUSIONS 

We have presented a detailed spectroscopic and photometric study 
of OGLE 05155332-6925581 including the analysis of high res- 
olution spectra and the application of a sophisticated light curve 
model. The multi-band light curves covering 13.7 years of photo- 
metric data from the MACHO and OGLE projects were analyzed 
along with the radial velocities of both stellar components of this 
eclipsing binary. 

We find that the system is best modeled with a geometrically 
thin and optically thick disc around the gainer. The analysis of the 
photometric data has allowed us to derive improved orbital param- 
eters and physical stellar properties such as stellar masses, radii, 
luminosities and the effective temperatures. All these parameters 
are given in Table 4. Interestingly, the orbital solution remains con- 
stant during the long cycle. This means that the disc is not the origin 
of the long cycle, which is also consistent with previous claims that 
the source of this variability is not eclipsed (M08). On the other 
hand, the presence of variable Balmer emission suggests the exis- 
tence of a fourth component in the system with physical conditions 
usually associated to optically thin circumstellar gas. Due to the 
presence of this disturbing component, our line profile investiga- 
tion is not conclusive regarding gainer critical rotation. 

Taking advantage of our new and improved system and 
stellar parameters, we explored the evolutionary stage of 
OGLE 05155332-6925581 with the aid of pubUshed grids of evo- 
lutionary routes for binary systems of similar masses, considering 
conservative and non-conservative evolution. In particular, the grid 
of models by Rensbergen et al. were considered in our work. We 
find the system in a semi-detached configuration and at a stage of 
rapid mass transfer. The donor has exhausted its hydrogen in its 
core and we observed the systems in a Case-B mass transfer, with 
an age of 4.76 x lO' years. Contrary to the case of V393 Scor- 
pii, the gainer parameters do not deviate too much from a main 
sequence star of similar mass. This is consistent with a younger ac- 
creting object, in the sense that OGLE 05155332-6925581 has not 
had time to accrete large amounts of matter since it is still inside 
the burst of mass transfer (Fig. 8). On the contrary, V393 Scorpii 
already passed this event in the Ml 2a model and had time to accu- 
mulate an important amount of mass in a massive disc. Our studies 
indicate that the DPV phenomenon is observed at different evolu- 
tionary stages; Case-A for V 393 Sco (M12a) and DQ Vel (Barria et 
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Figure 8. Left. M2 (upper curve) and M2 for the best evolutionary model, X'^ is shown in the inset graph. The few spikes in the M2 curve reflect minor 
convergence artifacts produced during the numerical calculations. Right. Evolution of orbital period and mass of the components with the mass transfer time 
for OGLE 05155332-6925581 with the initial orbital period Porb.i = 2.5 days. The vertical dashed lines and filled circles indicates the position for the best 
model. 



al., in preparation) and Case-B for OGLE 05155332-6925581 but 
always inside or after a main mass transfer burst. 

One notable observation in OGLE 05155332-6925581 is the 
absence of orbital period change, overall considering the relatively 
large mass transfer rate. This cannot be explained in terms of the 
mechanism of mass and angular momentum loss proposed by Rens- 
bergen et al. (2008, 2011). Actually, the best model indicates that 
the system should be found at a phase of fast period increase. We 
argue that mass outflows through the Lagrange L2 and L3 points 
(as claimed by MOB to explain infrared spectroscopic observations) 
could extract angular momentum from the system in order to bal- 
ance the mass exchange and produce a relatively constant orbital 
period. In addition, if the Balmer and He I emissions are due to a 
bipolar wind as suggested in the case of the DPV V393 Scorpii 
(Ml 2b), this wind could be an additional channel for escape of an- 
gular momentum from the system. 

It is worthy to mention that current theoretical evolution- 
ary models for close binaries do not include neither predict 
the DPV phenomenon. As the number of DPVs is relatively 
large, the phenomenon corresponds to a relatively long phase 
in the life of close binaries of intermediate mass. Our study of 
OGLE 05155332-6925581 suggests that the DPV phenomenon 
could have an important effect in the balance of mass and angu- 
lar momentum during the system evolution. This makes sense if 
the long cycle in OGLE 05155332-6925581 turns out to be due to 
a cyclic bipolar wind as suggested for the DPV V 393 Scorpii by 
Ml 2b. 
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